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■ Abstract Stationary differential systems witli polynomial right sides are considered. 



Necessary and sufficient conditions are formulated when a given domain is a domain of 
asymptotic stability and the origin of coordinates is either focus or center. The problem 
of construction of a stabilizing control in the form of polynomial is studied. 

Key words Differential systems of n-th order, asymptotic stable systems, stabilizing 
control, domains of asymptotic stability. 
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In the lO*'^ century the great French mathematician Henri Poincare formulated the problem 
^ ! of finding stability conditions for differential systems without calculation being a solution. 

^ I At the end of 19*^ and at the beginning of 20*^* centuries the great Russian mathematician 

A.M.Lyapunov developed the mathematical stabilization theory for differential systems. 
For this purpose he developed two methods. The first one is based on the characteristic 
numbers of the fundamental matrix. The second one is based on construction of special 
functions called Lyapunov's functions which have properties similar to distance from a 
considered point to the origin of coordinates (it is supposed that the zero solution is 
stationary). Lyapunov's methods received world recognition. 

In this paper, differential systems with a right polynomial part /(.) of 

X = {xi,X2, e 

are considered i.e. 

/(X) = (/1(X),/2(X),...,/„(X))*, 

where 

JpK-i-J — 2^ "/i,/2,...,in 1 2 ■••■^n ) 
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(p) 

p G 1 : n,li,l2, ■■■In arc non-negative integers and * is the transposition sign, al^i^ are 
real- valued numbers, Ip is the set of degrees of the polynomial fp{x). 
To analyze the system 

x= f{x) 

a method is suggested which is different from Lyapunov's methods and based on a system 
transformation idea, so that we arc able to say something definite about stability. 

From the technical point of view the important problem is finding a domain of asymp- 
totic stability. Conditions are given so that any domain including the origin of coordinates 
is one of them. Necessary and sufficient conditions are found when the origin of coordinates 
is either focus or center. 

iv) 

The case is considered when the coefficients o, , , depend on t. Sufficient conditions 
are given for the problems formulated above. 

The systems considered are interesting because asymptotic stability of the zero solution 
is equivalent to the following statement, that any solution of our system starting from any 
point of some small region of the origin of coordinates tended to be zero; this is not correct 
in general case. 

Further, we solve the problem of constructing a stabilizing control in any given region 
of the origin of coordinates for the system 

x= f{x,u), 

where x £ W" isa phase vector ,u & R^' isa control, f{x, u) = {fi{x, u), f2{x, u), fn{x, u))* 
is a vector polynomial of x and u i.e. 



fp{x,u)= 

h,l2,---,ln,mi,m2,...,mr€lp 



"■h,h,...,ln,mi,m2,...,mr-'''l-^'2 ■■■■''n H "2 ••■"r- 



(p) 

p £ 1 : n, /i,Z2,...,^n,"T.i, 771-2,..., are non-negative integers, a.i\^i2^..._/„^„j^„2^...^„^ are 
real- valued numbers, Ip is the set of degrees of the polynomial fp{x,u). We assume that 
the zero vector = (0, 0, 0) G 3fi"' is a solution of the system i.e. /(0, 0) = 0. 

The theorem was proved that any region including the origin of coordinates can be 
made a region of asymptotic stability if we choose a suitable control u{x). The control u{x) 
is chosen as a polynomial with a degree not higher than the degree of the vector-function 
f{x,u) as a function of x. 

The formulated problem about asymptotic stability of a system with a right polynomial 
part is very important for practical applications in physics and technics and can be found 
in the book by V.I.Zubov "The lectures on the control theory", 1975. p. 60. 
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2 Domains of asymptotic stability 



Let us consider the differential system 

^= fix), (1) 
where x = {xi,X2, ■■■,Xn) £ and 

f{x) = {h{x),f2ix),...,fn{x)r. 

The vector-polynomials fp{.) and their coefficients af^ 

,h,-,ln ^^^i^fy t'^^ conditions in the 

Introduction. 

We assume that /(x) / for all x 7^ in some neighborhood of the origin of coordinates 

0. 

We will call the order deg{fp{.)) of the polynomial /p(-) the maximal degree of the 
polynomial fp{.) in the variables Xj,j G 1 : n, in totality i.e. 

deg{fp{x)) = max {h + I2 + ... + In)- 

So if for n = 2 and lij 7^ 

fi{x) = luxj + li2xl + I1ZX1X2, f2{x) = hixl + I22X2 + I23X3, 

then the order of fi{x) is equal to two and the order of f2{x) is equal to three. 

We will call the order of the function /(.) the maximal degree of the polynomials 
fp{x),p G 1 : n, regarding the variables Xj,j G 1 : n, i.e. 

degifix)) = max deg{fp{x)). 

pel:n 

Consequently, the order of the function /(.) for the example written above is equal to 
max(2, 3) = 3. 

We assume that the definitions of stability and asymptotic stability are known (see 
[l]-[4]). 

Definition 2.1 [1], [2] A domain D,0 ^ intD, consisting of whole trajectories of the 
system (1), is called a region of asymptotic stability if the limit 

\\x{t,XQ,to)\\^Q (2) 

when t 00 is fulfilled for any initial point xq E D and any solution x{-), x{to) = xq, of 
the system (1). In this case we will call the system (1) asymptotic stable in D. 

The following problem is important for practical applications. 
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Problem 2.1 It is known that zero solution of the system 



(3) 



A[n xn\,x E is asymptotic stable if all eigen-values of the matrix A have negative real 
parts. In this case the system (3) is asymptotic stable in 3?". 
Consider now the differential system 



where ip{-) is a vector polynomial with degree not less than two, = 0. 

What are the conditions on the coefficients of the vector- polynomial ip{-) from (4) 

under which the system (4) is asymptotic stable in a given domain D,0 £ intD, i.e. for 
any solution of the system (4) and any initial point xq E D the limit (2) is true? 

The considered circle of questions includes the problem about center and focus. This 
problem is formulated in the following way. 

Problem 2.2 Assume that the origin of coordinates (0,0) for the system (3). n = 2, 
x € 3?^ is the center. The necessary and sufficient conditions for this are that all eigen- 
values of the matrix A are imaginary. By adding a vector polynomial ip{-) (system (4)) 
the center (0,0) can be a focus. It is needed to find conditions on the coefficients of the 
polynomial ip{-) that the point (0,0) was the center of the differential system (4). 

Let us rewrite the system (1) in equivalent form 



The elements aij{x) of a matrix A{x)[n x n] are continuous polynomial functions of x. 
Conversion (4) to (5) is not unique. It can be done in an infinite number of ways. In fact, 



x= Ax + ip{x) 



(4) 



x= A{x)x. 



(5) 



aij{x) = 



E 



a 




{x)x['x^^...x] 



...X, 



n 



li+l2+-+ln<deg{f) 



if Ij > 1. We have the following correlation for the coefficients 




(6) 



for all X from some region D,0 E intD. 

For instance, let the system (1) have the right part 




Then 




aiXiX2 + 02^2 1 + PlXi + P2X1X2 

■ J1X1X2 + J2X2 SiXi + S2X1X2 
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where the coefficients ai, f3i,^i, f5i,i = 1,2 such that 

ai + A = 1, 71 + = 3, 

a2 + P2 = 2, 72 + (^2 = -2. 
The system (4) for the given vector function /(•) can be rewritten in the following form 

x= Aqx + C{x)x (7) 

where 

.2 a ^2 



\ _ I aiXiX2 + 023^2 l^l^i + 1^2X1X2 
V 713^1^2 + 722^2 5ixl+ 82X1X2 




Since the eigen-values of the matrix Aq are Ai^2 = the point (0,0) is the center for the 
linearized system. But it is question for the nonlinear system (7). 

We will consider all possible continuous matrices A{x) whose elements arc polynomial 
functions of x for which the system (5) is equivalent to the system (1) . We will denote 
the set of all such matrices by A. 

Let us solve the problem 2.1. 

Theorem 2.1 In order that a domain D consisting of the whole trajectories of the system 
(1) i.e. x(-,xo,io) G -D,xo G -D, for all t > to is a region of asymptotic stability it is 
necessary and sufficient that there is a matrix A{-) E A of the system (5) in the domain 
D whose eigen-values have negative real parts at any point x E D,x ^ 0. 

Proof. Necessity. Let a domain D be asymptotic stable. Consider any trajectory 
x{-,xo,tQ),x{to) = xq. There is such a transformation ^ = X{x) of in a neighborhood 
of a point xi = x{ti) that transforms the system (1) to the differential system 

e= B{o^ (8) 

where the matrix B has the eigen-values with negative real parts. 

Indeed for any point x{ti,xo, to),ti > to there is a linear transformation defined in a 
small neighborhood of the point x{ti,xo,to) and that is the linear part of the transforma- 
tion X{.) in this neighborhood, so that the system (1) is transformed by this transformation 
to the diff'erential system (8) and the vector 

Xi =x {ti,Xo,to) 

is transformed to the vector 

6 =C {ti,xo,xo). 

Any differential equation can be defined by its current of tangent vectors. If we choose an 
asymptotic stable linear differential system like the system (3) which current of tangent 
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vectors is close to current of tangent vectors of our system locally for each point x of the 
domain D then the eigen-values of the matrix of this asymptotic stable linear differential 
system will have negative real parts. Consequently, the above mentioned transformation 

X(x) exists. 

As soon as the system (8) is transformed to the system 

x= T~^B{^i)T^x 

under the linear transformation 

in a neighborhood of the point and all locally linear transformations of the system (1) 
have the form (5), then the matrix A{-) G A exists for which 

A{x{h,xo)) = T-'B{^i)T, (9) 

and the eigen-values of the matrix A{x{ti)) are equal to the eigen-values of the matrix B. 

Since the transformations and T"^ are continuous with respect to x, the matrix 
A{.) is also continuous with respect to x £ D . 

Sufficiency. Let the system (1) admit such a transformation to the form (5) that all 
eigen-values of a matrix A{x) have negative real parts for all x G D, x / 0. We prove that 
the domain D is a region of asymptotic stability. 

Fist we remark that in general case not any solution of the system (1) 

x{t) = x{t, Xq, to), x{to) = Xo 

can be represented in the form 

• , . r* A{x{T))dT , . 

X [t) = e-'*o Xq (10) 

where integration is taken along the integral curve x{.). 

Instead of (10) we will use the sequential approximations {xk{t)} of the system (1) on 
the segment [to,t] having the form 

a;fc(t) = e^(^'=-i)(*-*'=-i)xfc_i(t) (11) 

where {ti} is a subdivision set of the segment [to,t] and 

Xi+i = e^(^^)(*^+i-**)a;i, i G : (/e - 1). (12) 

If we substitute (12) into (11) we will get 

xi^(t) = e^(^'^-i)(*~*'^-i)e'^^^'^-2^'^*'^-i~**'-2^...e'^^^°^'^*i~*°^xo. (13) 

Since A{xi) — > ^ _^ ^^^^ A{x{t)), where x{.) is a solution of (5), we have 
ti T 



Xk{r) =^ x{t) 



(14) 
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uniformly on r G [to, t] when — >■ oo and 

max I ti — ti-i |— >-fc 0. 

It is obvious that all matrices A{xi),i G : (A; — 1), have the eigen-values with negative 
real- valued parts for enough big k. 

Denote by Xi{xj) the eigen-values of the matrix A{xj),i & 1 : n,j & : k — 1. Then a 
number C > exists that 

||gA(a;fe_i)(t-tfe_i)gA(a;fe_2)(tfe_i-tfe_2)_ _gA(a;o)(ti-to)^p|| < C'e'^^*'"^*~*°-' ||xo || 

where 

X{t) = lim max Re Xi{xj). 
fc->oo j G : A; - 1, 

i E 1 : n 

Since X{t) < for any t > to, 

X{t){t-to) < 0. 

Let 

-a < X{t){t-to) < 0,a > 0. 

This relation means that the trajectory x{-) can not go to a stationary point or be a 
stationary orbit. Indeed if the trajectory x{.) has a stationary point or is a stationary 
orbit then 

lim X(t)(t — to) = — oo. 

t— >oo 

Consequently, 

II x{t,xo,to) II — > 

when t — > oo. We obtain the contradiction. Thus, the origin of coordinates can be only 
a stationary point. The sufficiency and the theorem are proved. □ 

Remark 2.1 It follows from the Necessity of the Theorem 2.1 that there is one-to- one 
correspondence between the matrices A{x) & A of the system (5) and locally linear trans- 
formations Tx of the system (1) defined in a neighborhood of the point x(see (9)) . Indeed, 
there is a matrix A{.) for any transformation and, on the contrary, there is own coor- 
dinate system for any matrix A{x) G A, in which the equation (9) is true. 

Corollary 2.1 . It follows from the Necessity of the Theorem 2.1 that the degree on x of 
any element of a matrix A{x) which we denote by aij{x),j G 1 : n, that is a polynomial of 
X, does not exceed the degree of the vector-function /(•). 

Proof The sum of the coefficients aij{x) ■ Xj over j is equal to zero for those elements 
that do not belong to /«(•). And on the contrary the sum of the coefficients ce^j]^ 
over j is not zero for those elements that are the terms of fi{ ). It follows from here that 
choosing a"^\^ from the system (6) we can only consider the elements of the vector- 

polynomial /(•). □ 
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Corollary 2.2 . For asymptotic stability of a solution x{.,xo,to), x{to) = xq, of the 
differential equation (1) it is sufficient that there was 6{to,e) > for any to,e > that for 
II ^0 ||< <^ 



Proof. As soon as the eigen-valucs of a matrix A(-) have negative real parts in any 
neighborhood of the origin of coordinates from which any solution tends to zero-vector , 
then for any k we have from (13) that 



i.e. all solutions are in an e -neighborhood of the origin of coordinates. The latter means 
stability. The Corollary is proved. □ 



3 To problem about center and focus 



The following problem is interesting from technical point of view: to recognize focus or 
center i.e. to give the conditions when trajectories turn around and go to the origin of 
coordinates or remain closed. 

From the start we will give the definitions of focus or center. Beforehand we define 
trajectories turning around some ray with the initial point G infinitely often. 

Definition 3.1 Let us say that a trajectory x{t,xo,to) turns around a ray I G 3?" with the 
initial point G infinitely often if the radius vector r{t) = x{t,xo,to) forms with the 
ray I an angle ip{t) taking all values from a segment [</'i,¥'2] infinite number of times. 

Definition 3.2 We will call the point G 3?" focus of the system (1) if 

1) the point is asymptotic stable ; 

2) any solution x(i,xo,io) of the system (1) turns around a ray Z G Ji"" with the initial 
point G 3?" infinite number of times. 

Definition 3.3 We will call the point a center of the system (1) if 

1 ) the point is stable ; 

2) all solutions x{t,XQ,tQ) of the system (1) remain in some neighborhood of the point 
and are closed loops. 





Xk{t) \\<\\ xo \\< e 



and in limit on k 



x{t,XQ,tQ) \\<\\ xq \\<e 



9 



Definition 3.4 We will call the point G 3^*^ a node of the system (1) if 

1 ) the point is asymptotic stable ; 

2) there are not a ray I G di'" with the initial point G and a solution x{t,xo,to) of 
the system (1) turning around the ray I infinite number of times. 

Theorem 3.1 In order that the point E is a focus of the system (1) with a right 
polynomial part /(.), f{x) ^ for x it is necessary and sufficient that 

1. There is such a continuous matrix A{-) G A of the system (5) that all its eigen-values 
at any point x,x ^ 0, from some neighborhood D of the origin of coordinates, con- 
sisting from the whole trajectories, have negative real parts and non-zero imaginary 
parts. 

2. There is no matrix A{-) £ A of the system (5) with the negative real-valued eigen- 
values at all points x & D,x 0. 

Proof. Necessity. Let the point be a focus for the system (1). There is such a 
transformation ^ = X{x) so that the system (1) can be rewritten in a neighborhood of 
some point xi,xi ^ 0, in the form 

where is a matrix whose eigen-values have negative real parts and nonzero imaginary 
parts. 

Consider any trajectory x(-, xq, to),x(to) = xq. Then for any point xi G x{t, xo,to),xi / 
0, there is a transformation T-cj in a neighborhood of the point Xi that the system (1) can 
be rewritten in the form indicated above. 

If the trajectories of the system (1) do not turn around any ray I G 3?" with the initial 
point G Sf?*^ infinite number of times then such system can be transformed by some 
continuous transformation that corresponds to some matrix A{-) G A that is the matrix 
of the system (5) whose eigen-values are negative real- valued numbers at all points x E D 
(theorem 2.1). The latter is impossible according to the Condition 2 of the theorem. 

Sufficiency. Let the conditions 1,2 of the theorem be true. Prove that the point is 
a focus of the system (1). 

Any solution x{-,XQ,tQ),x{to) = xq, is the limit of some sequence x^.(t) on /c — > oo 
obtained from (13). According to the conditions all matrices A{xj),j G 1 : A;, have the 
eigen-values \i{xj),j G 1 : k,l G 1 : n with negative real- valued parts ai{-) and nonzero 
imaginary parts bi{-) : Xi{x) = ai{x) -\- ibi{x),i'^ = —1, Z G 1 : n, that are nonzero at any 
point x G -D. It follows from here that the point is either focus or node. 

We will prove that under the Condition 2 all trajectories turn around a ray I G 
with the initial point G ^R."" infinite number of times. If it is not true, then there is such 
a transformation which corresponds to a matrix A{-) G A whose eigen-values are negative 
real-vahicd numbers ior x E D,x ^ 0. This contradicts the Condition 2. The sufficiency 
and the theorem are proved. □ 
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Corollary 3.1 If instead of the Condition 2 of the Theorem, 3.1 we require that there is a 
matrix A{.) E A of the system (5) with the eigen-values Aj(x), i £ 1 : n, described in the 
condition 1 for which 

Im \i{x) > a, a > 

or 

Im Xi{x) < a, a < 0, 

where Im\k{.) denotes the imaginary part of X^i-), for x £ D,x ^ 0, then all trajectories 
turn around a ray I G with the initial point € 3?"" infinite number of times. 

Proof. From the representation of any solution of the system (1) as a limit of sequence 
{xj{.)} from (13) it follows that a solution x(-,xo,to) has a finite number of turns around 
any ray / G SR" with the initial point G 3?" if and only if for any i e 1 : n the sum 

Yl IrnXi{xj-i){tj -tj-i), 

jel:k 

where Xj{.) was defined in (13), or in the limit on ^ oo the integral 

roo 

/ Im Xi{x{T))dT 
Jto 

is convergent. But under the condition of Corollary 3.1 this integral is divergent. The 
corollary is proved. □ 

Mathematicians were concerned about the problem of recognition of center ore focus for 
a long time. The theorem given below for two dimensional system with a right polynomial 
part has the point = (0, 0) as a center. 

Theorem 3.2 In order that the point = (0, 0) be a center for the two-dimensional 
system (1) it is necessary and sufficient that there was a matrix A{.) G A of the system 
(5) whose eigen-values are non-zero imaginary numbers in a neighborhood S of the origin 
of coordinates, consisting from whole trajectories, where f{x) ^ for x ^ Q,x E S. 

Proof. Necessity. Let the point = (0, 0) be a center of the system (1). Prove that the 
condition of the theorem is true. 

Take such a small neighborhood S of the origin of coordinates where all trajectories 
starting in S are closed loops. Consider any trajectory x{-, xo,to), x{to) = xq. Take a point 

Xi = x(ti,XQ,to). 

There is such a transformation in a small enough neighborhood of the point Xi,xi ^ 0, 
that the system (1) can be rewritten in the form 

^= Bm 

and the matrix B{.) has the imaginary eigen-values for all ^ G 5". As was proved above 
(Theorem 2.1, the Proof of the Necessity) the matrix B{-) corresponds to some matrix 
A{-) that like the matrix B{-) has the imaginary eigen-values. The necessity is proved. 
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Sufficiency. Let the condition of the theorem be true. Prove that the point = (0, 0) 
is the center. 

As soon as any transformation Tx is defined in a small neighborhood of the point x, 
then the above statement is not sufficient for the point = (0, 0) to be a center. The point 
= (0, 0) may happen to be a focus convergent (nonconvergent) with a finite number of 
turns around any ray I G 3f?" with the initial point 0. But it is impossible because from 
the limit of the sequence (13) we conclude that for any k 

II Xk{t) 11 = 11 Xo II . 

From (14) it follows that the solution x(-,xo,to) can't go to the origin of coordinates. 

A set of closed and not closed loops could occur. From the representation of an integral 
curve x{-,xo,to) in the form (13) and (14) it follows that there is not a matrix A{-) G A 
with the imaginary eigen-values at all points of a small neighborhood of the origin of 
coordinates. The latter contradicts the condition of the theorem. 

There is another case when for all i G 1 : n the sum 

^ ImXi{xj-i){tj - tj-i) 

jel:k 

or in the limit on A; — > oo the integral 

roc 

/ Im \i{x{T))dT 
Jto 

has a finite value where Aj(-), i e 1 : n, are the eigen-values of the matrix A{-). In this 
case the trajectory x{-,XQ,tQ) goes to some stationary point x G S*, a; 7^ 0. It means that 
f{x) = 0. This contradicts the condition of the theorem. Consequently, the point is the 
center. The sufficiency and the theorem are proved. □ 

Corollary 3.2 // a matrix A{.) £ A of the system (5) exists with the imaginary eigen- 
values Xi{x) for all X & S where S is a neighborhood of the origin of coordinates 0, 
consisting from whole trajectories, and 

Im Xi{x) > 

or 

Im Xi{x) < 

for all X £ S,x ^ 0, and i E 1 : n then the point is the center of the n-dimensional 
system (1). 

Proof. Indeed, if the conditions of the Corollary 3.2 hold, then the conditions of the 
Theorem 3.2 hold as well. □ 

The Corollary 3.2 does not require that f{x) 7^ for all x 7^ from some neighborhood 
of the origin of coordinates. 
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Remark 3.1 The theorem can be proved for any n- even- dimensional spaces. It is not 
difficult to see that for n dimensional spaces with n odd, n = 2k-\-l, k is a natural number, 
the system (1) can not have the origin of coordinates as a center. Indeed, there is no 
imaginary number among the eigen-values of any matrix A(-) for odd n. 

Remark 3.2 Let us consider the differential system 

x=f{x,t) (15) 
where x = {xi,X2, Xn) G 3ft" and 

f{x,t) = {f^{x,t),f2{x,t),...,fn{x,t)y. 

The vector-function fp{.) is a polynomial of x i.e. 

ll,l2,---,ln&Ip 

p G 1 : n,li,l2, ...In are non-negative integers, a\^\^ i^it) are continuous real-valued func- 
tions and * is the transposition sign, Ip is a finite set of indexes of the polynomial fp{.). We 
will assume that f{x, t) for all x ^ in some neighborhood of the origin of coordinates 
and t > to. 

In this case we will denote by Xi{x,t),i £ 1 : n, the eigen-values of a matrix A{.). If 
we demand that all cited above statements about the eigen-values Xi{x,t),i G 1 : n, are 
true for all t > to then we obtain the sufficient conditions for the cited above theorems and 
corollaries. 

Problem 3.1 For the system 

x=f{x,t) (16) 

where /(., .) is a polynomial of x and t it is required to find some conditions when the 
system (16) is asymptotic stable in a domain D,0 £ intD. 

The ideas stated above do not apply to the system (16) because the terms of /(., .) 
can be unbounded along some solution x{-,xo,to) in D. If we require the terms of /(.,.) 
to be bounded along the trajectories of the system (16) and that {t,xo,to) \\ goes to 
zero uniformly on xo when t — > oo, then this system can be transformed to an equivalent 
stationary system. This idea will be developed in following articles. 



4 Stabilizing control 



Let us go to the question of finding a stabilizing control. 
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Let us consider the differential system 

x=f{x,u) (17) 

where x = (xi, X2, x^) € 3(?" is a phase vector, u = {ui,U2, ■■■,Ur) E W is a control, 
f{x,u) = {fi{x,u), f2{x,u), fn{x,u))* IS a vector-polinimial of x and u with constant 
real-valued coefficients, i.e. 

fp{x,u) = ^ (^h\2,---,lri,mi,m2,...,mr^l ^2 ■■■^n^T^ 

(p) 

p e 1 : n, Ii,l2,...,ln,mi,m2,...,mr are non-negative integers and a;\,i2,...,/„,^j,^2,.-,'nr 
are real- valued numbers, Ip is a finite set of indexes of the polynomial fp{-)- 

Let us assume that the zero- vector = (0,0,..., 0) G 3^^"' is a solution of the system 
(17) for « = € 3fi^ 

Definition 4.1 The control u{x) = {ui{x),U2{x), ...,Ur{x)) G K'' is called stabilizing 
in a domain D £ 3?", G intD, for the system (17) if any solution of (17) x{t) = 
x{t,xo,to,u{x{t))) satisfies the limit 

II x{t,XQ,to,u) II — > (18) 

when t — y oo, x(to) = xq G D. 

As stated above the condition (18) is sufficient for the system (17) to be asymptotic stable 
in D, i.e. the zero solution of the system (17) is asymptotic stable and the limit (18) is 
true for any initial point x(to) = xq E D. 

Problem 4.1 It is required to find such a stabilizing control u = n(x) in a given domain 
-D, G intD,x G D that for any solution x{t,xo,tQ,u),x{to) = xq of the system (17) the 
correlation (18) was true. 

For linear system 

x= Ax + Bu, (19) 

where B[n x r] is a matrix of amplification coefficients, a stabilizing control w G 3fi^ can 
be chosen in the form u = Cx so that (18) is true. 

Theorem 4.1 [lj,[4j- If the rank of the system 

B, AB, A'^B, A^'-^B 

is equal to n then we can always construct a stabilizing control in 3^?" in the form 

u = Cx, 

where C[r x n] is some matrix. 
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Since the matrix B is defined by technical essence of the system, we can choose it con- 
structing a system by ourselves . Therefore the condition of the theorem can be always 
satisfied. 

The theorem was written for analogy and comparison with the results below. 

Our goal is to construct a stabilizing control u{x) that a domain D were a domain of 
asymptotic stability. In common case it is not always can be done. But if we change the 
system (17) a little bit the problem can be solved. Instead of the equation (17) we will 
consider the equation 

x= f{x,u) + ^{u) (20) 

where ip{-) is a vector-polynomial (fiu) = {ipi{u),ip2{u), ...,ipn{u))* with the degree not 
bigger than the degree of the function f{z) as a function oi z = {x, u), 

(p) 

where p E 1 : n and b^^^^ are constant real- valued numbers, zi, Z2, v are non-negative 

integers, Mp is a finite set of indexes of the polynomial V'p(-)- 

In practice it is possible to construct Lp{-) because wc choose a control n(-) ourselves. 
We will look for a stabilizing control u = u{x) in form of polynomial in x. 

Theorem 4.2 For any domain D, G intD, the vector-polynomials u{.) and (p{.) can be 
chosen such that 

1) D is a region of asymptotic stability for the differential system (20); 

2) the degree of u{x) does not exceed the degree of the vector-polynomial f{x, u) as a 

function of x; 

3) the degree of ip{.) is not bigger than the degree /(.) as a function of 
z = (x, u). 

Proof. We will prove that a vector polynomial u{x) can always be chosen satisfying the 
following conditions 

1) the degree of u{x) does not exceed the degree of the vector-polynomial f{x,u) as a 
function of x; 

2) u(x) is a stabilizing control for (20) in the domain D. 
We will use the results obtained before (Theorem 2.1). 

Let us substitute in (20) the control u = u{x) in the form of a vector-function of x. 
The system (20) is rewritten as 

x= fix) (21) 

or 

x= A{x)x. 

We can write conditions for the matrix ^4(0;) to have the eigen-values with negative real 
parts. 

There is no difficulty calculating the degree of the function f{x,u) and the function 
(p{u) as a function of x after substituting u = u{x). 
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Let us denote by Ix and lu the degrees of the function f{x,u) in the variables x and u 
correspondently. Then after substituting u = u{x) the degree of the function /(.,.) as a 
function of x is not bigger than Ix + lulx and the degree of the function (p{.) as a function 
of X is not bigger than lx{lx + lu)- It is easy to see that 

Having chosen the function '^{■), the coefficients of the vector polynomial f{x) can be 
chosen so that the matrix A{x) has the eigen-values with negative real parts. The theorem 
is proved. □ 

In an analogous way it can be proved the following theorem. 

Theorem 4.3 For the differential system (20) with even n the vector-polynomials u{-) 
and (/?(•) can he chosen such that 

1) the origin of coordinates is the center for the system (20); 

2) the degree of u{x) does not exceed the degree f{x,u) as a function of x 

3) the degree of (p{.) is not bigger than the degree of f{.) as a function of z = {x,u). 

This theorem can be used in physics of plasma specially for stabilization of plasma in 
reactor. 

(v) 

Remark 4.1 In the case when the coefficients a^^^^ mr depend on t the coef- 

(p) ! ) 1 1 ' > > 

ficients of^ ^^ will depend on t as well. We can try to find these coefficients that the 
conditions of the theorems of the previous section were true for all t > to. 



5 One aspect of application 



We will study for instance how to choose a control u(-) so that a given domain D consisting 
from the whole trajectories of a differential system were a region of asymptotic stability. 
For that we have to solve the following optimization problem. 

Let us substitute a vector-control u{-) in a form of polynomial of x with a degree m 
not bigger than the degree of the the vector-polynomial f{x,u) as a function of x in the 
equation (20) and rewrite (20) in the form (21). Denote the eigen-values of a matrix 
i(-) in (21) by Aj(x,aij-,4f,j2,...,j„,?'lf,i2,...,iJ where and 4I!i2,...i„ ^re the coefficients 
of the matrix A{-) and the vector-function u{-) correspondingly, b\^\^ ^ , p G 1 : ra, are 
the coefficients of the vector- function </'(.). Then our problem is reduced to the following 
optimization problem: to find out such continuous functions aij{x),i,j £ 1 : n, x E D and 

the numbers , b^\„ , p G 1 : n, that the following correlation 

ReXjix,aijix),d$^,^^___^,^,b^,^_,J < (22) 
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was true for all x & D,x 0, where aij{x),i,j G 1 : n, are connected with each other by 
hnear equations for each x e D,x 0, hy other words we should solve the next problem 

sup inf Re Xj{x, aij,dfli ^^,1^^ J < V (5 > 0, 

where -B^(O) = {z £ 3?" 1 1| z \\< S}, S— is any sufficient small number for which D D i?^(0). 

It is obvious that the inequalities (22) can be replaced by an equivalent system of 
inequalities for coefficients of the characteristic polynomial of the matrix A{.). To solve 
this system it is easier than to find the eigen- values of the marix A{.). 
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Abstract In the article the necessary and sufficient conditions for a representation of 
I Lipschitz function of two variables as a difference of two convex functions are formulated. 

' An algorithm of this representation is given. The outcome of this algorithm is a sequence 

of pairs of convex functions that converge uniformly to a pair of convex functions if the 
conditions of the formulated theorems are satisfied. A geometric interpretation is also 
given. 
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1 Introduction 



The solution to this problem is very interesting for the fields of optimization ( for quasid- 
ifferential calculus [1]) and geometry. In [2] the author started to investigate the surfaces 
that are graphs of functions representing as a difference of two convex functions (so called 
RDC functions) to construct the inner geometry of such surfaces. This problem has been 
investigated by many authors (for example see [3] - [7]). 
■ The necessary and sufficient conditions for such a representation of a function of one 

variable are well known. These conditions can be written in the following way. 

Let be X — )• f{x) : [a, 6] — )■ K a Lipschitz function. Define a set, where the function /(•) 
is differentiable, as Nj. Since the function /(•) is Lipschitz, Nj is the set of full measure 
on [a,b]. 

For the function /(•) to be represented as a difference of two convex functions it is 
necessary and sufficient that the following condition was true 

sup{^^}cNf I •^'(^«) ~ \= V(/';o,&) < oo 

i 

where the derivatives are taken where they exist. The sign V means the variation of the 
function /' on the segment [a,b]. 

In the article [7] the necessary and sufficient conditions for a representation of any 
one degree homogeneous Lipschitz function of three variables as a difference of two convex 
functions are given. This result can be extended to the m*^ degree homogeneous functions. 
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I will now ignore the condition regarding to homogeneity and will consider any bivariate 
Lipschitz function /(•) with a Lipschitz constant L: {x,y) — ^ f{x,y) : D — >■ K where D 
is an open bounded convex set in 3?^ so that its closure D is compact. I will describe a 
representation with an algorithm giving a sequence of pairs of convex functions. These 
convex functions converge uniformly on D to a pair of convex functions if the conditions 
of the theorems formulated below are satisfied. 

Let p{D) be a class of curves on the plane XOY such that bound closed convex sets in 
the region D. Any curve r G p{D) will be written in the natural parametric representation 
i.e. the parameter r of a point M on the curve r(-) is equal to the length of the curve r(-) 
between an initial point and the point M. I will denote such a curve hyr{t),t G [0, Tr] . With 
the help of the curves r G p{D) the necessary and sufficient conditions for a representation 
of the function /(•) as a difference of two convex functions can be written in the following 
way. 

Theorem 1.1 In order that a Lipschitz function z — f{z) : D — > K was represented as a 
difference of two convex functions it is necessary and sufficient that 

iMD, /) > 0)(Vr G p{D)) V 0, T,) < c(L>, /) (1) 

where = /(r(t)) Vi G [0,r^]. 

The algorithm of a representation of the function /(•) as a difference of convex functions 
is given and it is proved that the sum of the obtained functions converges if the condition 
of the theorem 1.1 is satisfied. 

The proof is based on a special algorithm of a representation of the function /(•) as a 
difference of two convex functions. As the result we will get finite or infinite number of 
convex functions, which we will add together. It follows from the algorithm that the sum 
of the variations of the derivatives of these functions along any segment [a, b] C D is finite, 
if (1) is satisfied, and, consequently, the sum of the convex functions, obtained from the 
algorithm, is bounded above. 

For a representation of /(•,•) as a difference of two convex functions I will use two 
operations, which give us an opportunity to produce finite or infinite number of functions 
with bounded connected supports. According to the definition, a continuous function with 
bounded connected support is a continuous function such that is not equal to zero only 
on an open bounded connected set Q C .D and equal to zero on the set D\Q. The closure 
of the set Q is called the support of this function. A function with such support is called 
a finite function. 

I will denote the set of all finite functions on D by S. 

The first operation is a construction of a maximum convex function that does not 
exceed the initial one at all points of the set D. Using this operation in iterative way and 
subtracting from the initial one the obtained convex function I will get finite functions with 
supports Di, Di C D. After obtaining a finite function at some step I will apply to it again 
the construction of a maximum convex function on the convex hull coDi of its support. 
After that I will subtract from the maximum convex function the finite function, for which 
this maximum convex function was constructed, and will obtain new finite functions with 
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supports Dij C coDi and so on. It will be proved that the obtained finite functions will 
trend to zero on the norm in C{D) uniformly regarding to the step of iteration. 

The operation of the construction of the maximum convex function that docs not 
exceed the initial function is the similar to the operation that was used by B.A.Zalgallcr 
in [5] . He investigated possibility of a representation of some special kind of a function as 
a difference of two convex functions. 

The second operation is an extraction of finite functions, obtained at the previous step, 
that are convex on their supports co Di, co Di C D,. My goal is to represent each of them 
as the sum of convex and concave functions defined on D. When the condition of the 
theorem 1.1 holds the sum of the obtained convex functions is bounded above by some 
constant on the set D. 

I will call functions, represented as a difference of convex functions, RDC functions, 
according to the terminology in [2]. 

2 The Proof of the Theorem 

We will start the Proof of the Theorem from a description of the Algorithm. 

2.1 The algorithm of a representation of a function as a difference be- 
tween two convex functions 

I will construct a maximum convex function {x,y) fi{x,y) : D ^ ^ satisfying the 
inequality fi{x,y) < f{x,y). The function •) is Lipschitz with the constant L as soon 
as all its generalized gradients [8] are the gradients of the function /(•, •) at some points 
where the graphs Tf and Fj^ of the functions /(•) and /i(-) respectively touch each other. 

Using the extraction operation of finite functions from the difference /i (•,•)"/(■'■)) 
we can construct the finite functions (x, y) — >■ /i^jj (x, y) : D — >■ 3? with their supports 
Di^i-^ C D, ii G Ii C N (see fig.l). The sets Di^i-^, ii G Ii, are not necessarily convex and 
can have the common points on the boundaries, and also possibly that Ii is an infinite 
set of indices. The sum of the functions •) on ii is •) — /(•, •). If the sets -Di,ii 

are not convex, then I will obtain at the following step the convex functions on co Di^i^ , 
i\ G /i, by using the operation of construction of maximum convex function. The sets 
coDi/i^, ii e Ii. can have the common inner points. It will be proved below that •) 
is the Lipschitz function with the constant 2L. 

Our goal is to represent as a difference of two convex functions on D. If 

/i,n(")') is convex on its supports coDi^i^ then it can be extended to D in the following 
way (see [1], p. 118 ) 

/i,n(^)=max y^Qf^.^(^^^^ (/i,n (^) + ^ - ^))> 
z e coDi^i^ 

where z = {x,y),z = (x,y),dfi^i^{z) is the subdifferential of the function /i,ii(-, •) at the 
point z, as it is shown in [1]. The function /i,ii(-, •) is the maximum of the linear functions, 
consequently, it is convex [8]. 
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We can see that the difference 

is the convex function that is equal to zero on coDi^j^. From here 

fi,hi^) = fi,hiz) - h,n{z) G D, 

what was required to prove. 

If /i.ti(-,-) is not convex on -Di,ii then I will construct a maximum convex function 
— )• f2,ii{x,y) : coDi^i^ — > i? on coDi^j^ satisfying the inequality 

f2,h{x,y) < fi,n{x,y) y{x,y) e coDi^i^. 

The function f2,i^ (■> ") Lipschitz with the constant 2L on D as all its generalized gradients 
are the generalized gradients of the function •). The function /2,ti(-, •) is not equal 

to zero on intD2,ji = intcoDi^j^. Therefore the function /2,ji(-,-) can be represented in 
the same way as shown previously with respect to the difference of two convex functions 
on D: 

f2,h{z) = f2,h{z) - f2,h{z) '^zeD,z = {x,y). 

The difference f2,ii{z) — fi,ii{z),z G -D, is a Lipschitz function with the constant 2L and 
it will be proved later. The function f2,iiiz) — can be represented as the sum of 

the functions z /2,n,i2(^) '■ D ^ R with the supports -D2,ii,i2)^2 G h, respectively. In 
the general case the set I2 is an infinite set. Clearly that -D2,ii,i2 C coDi^j^. 
I will continue the described process further in the following way. 

The function {x,y) fk,ii,i2,...ikix,y) ■ D ^ R , that has been constructed at the 
k-th step has supports -DA;,ii,i2,-ifc ^ Dk-i,ii,i2,----ik-i- This function is defined on D and 
Lipschitz with the constant 2L. The set Dk,i^i2...ik C -D is its supports i.e. the closure of 
the bounded open connected set of the points where fk,h,i2,...ik{'^ ') is not equal to zero. 
The function fk,ii,i2,...iki'j') is equal to zero at other points that are not elements of the 
set Dk^ij^,i2,...ik- AH sets -Dfe,ii,i2, obtained during the described above process. It is 
possible that there is an infinite number of these sets. 

Description of the algorithm 

1. If at the k+1 - th step for some 11,12, ■■■ik the function fk,ii,i2...ik ('■> ') is not convex on 
its support D^^i^i^^^^^i^, I construct a maximal convex function {x,y) — >■ fk+iMM-'-iki^iV) '■ 
coDjfc^jj^j2,. -,«fc ~^ ^ satisfying the inequality 

fk+l ,ii,i2...ik 

{x, y) < fk 

The function fk+i,ii,i2,...ik{', ') is equal to zero at other points that are not elements of the 
set CO D^i^i^i^. 
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2. Let us define the difference 

I will represent it as the sum of finite functions, which I denoted by (x, y) — > fk+i,ii,i2,...,ik+i v) '■ 
D ^ R with bounded supports denoted by -Dfc+i,ii,j2,...,ifc+i C coDk/ti,i2,...,if,- 

3. I then let k=k+l and repeat the process again. In the general case the described 
process can give an infinite sequence of the finite functions fk,h,i2,---,iki'-> ')-> k & N. Accord- 
ing to the algorithm I have 

/(•,-) = /i(-,-)-E/Mi(-,-), 

h 

/l,n(') ■) = /2,n ~ /2,n,i2(') 

/2,ii,i2(') ■) /3,n,i2(')') " ^ ] /s, 11,12, 23 ')• 

h,i2,i3 

Continuing the expansion process we will have 
t 

/(•,•) = + X]^"-*-)*" X! /fc+l,n,i2,...,ifc(') •) + (-1)*"^^ X! /t+l,n,j2,...,it+i(-) •)• 
k=l 11,12, ---ik n,i2,---,it+i 

(2) 

I will show that all functions fk,ii,i2,...,ik{'i') Lipschitz with the constant 2L. For 
that it is sufficient to show that all generalized gradients of the functions fk,h,i2,—,ik 
are limited on the norm by the constant 2L. 

I begin the proof with one-dimensional case. 

At the first step after subtraction of the maximum convex function /i(-) from the 
function /(•) the gradients of an obtained function are equal to gi — f{x). The gradient 
g2 of the linear function (see the fig.2) becomes equal to gi — g2- It is obvious that 

lUi- V/(x) ||<2L, 

as gi is a generalized gradient of the function /(•) and || gi \\< L. 
At the second step we get a finite function with the gradients 

(5i-52)-(5i-V/(x)) = V/(x)-52 

in the region D2 and with the norm not bigger than 2L as far as || g2 ||< L. After two 
operations the gradient of a linear function in the region is equal to 

(51 - 92) - {g\ - 93) = 93 - 92- 
At the third step the gradients of an obtained finite function are equal to 
(53 - 92) - {91 - 92) + {91 - V/(x)) = 53 - V/(x), 
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in the region that is not bigger than 2L on the norm. 

We can continue our process further. Consequently, the gradients of all finite functions 
at each step are not bigger than 2L on the norm. 

In two-dimensional case the functions fk+i,ii,i2,...,i,,{'i ") can not be linear in the regions 
Dk+i,iii2....ik+i- We can always approach the function /(•) with the help of a piecewise 
linear function with any given precision and prove that all generalized gradients of the 
piecewise linear functions fk,ii,i2,...,iki'i ') limited on the norm by the constant 2L. 

Since any Lipschitz function /(•) can be approached by piecewise function with gra- 
dients of linear parts not bigger on the norm than L + (i where /x is any small positive 
number, it is necessary to follow how the gradients of linear functions change. We can 
see from computations, given above, that the gradients of the functions fk,h,i2,...,iki'^ ') 
are equal to zizg^ T where are the gradients of linear functions in the regions 

iii2----ik+i- follows from here that the gradients of the functions fk,ii,i2,...,iki'^ ') 
not bigger 2L on the norm. It is not difficult to prove it by the mathematical induction. 

I will prove that our process of the construction of the maximum convex functions that 

do not exceed the initial functions on their supports Dk.i1.i2 ik 1 ^'"^ '^ot stop while there 

are some points a, 6 G -Dfe,ii,i2,...,jfc with the different values of the function fk,ii,i2,...,iki'^ ') 
at the step k. Our process strives to make equal all values of the function fk,h,i2,-,iki'^ ') 
at all points of its support. 

This statement can be proved by the contraposition method. 

Consider at each step k the biggest deviation of the function fk,ii,i2,...,iki'j ') from zero. 
I will denote this by hk and the points where this deviation occurs by (xfc, yk) E D. Since 
D is the compact set, the sequence {{xk,yk)} has a limit {x, y) e D.I have already proved 
that the functions fk,ii,i2,...,iki'^ ') Lipschitz with the constant 2L. There are the points 
from D in any small neighborhood of the point (x, y) where the values of the functions 
/fc,ii,i2,...,ifc(") ■) s-i^s equal to zero and, consequently, differ from each other on the constant 
\ hk \> h > that contradicts the hypothesis of fk,h,i2,--iki.' ■> ') ^® Lipschitz with the 
constant 2L. 

Consequently, when k ^ 00 



uniformly on k. 

Moreover, all functions fk,h,i2,...,ik go to zero uniformly on k. This statement can be 
proved by contraposition in analogous manner like above. In the opposite case we would 
get a contradiction with uniform Lipschitz quality what can be seen from the previous 
considerations. 

Prom the algorithm and (2) it follows that the sum 



fk,H,i2,-,iki^^y) V(x,y) G D 





,1'k+l 




ll,l2,...,lk+l 



of the finite functions received at the k-)-l-th step also tends to zero when /c — )■ 00. For 
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any e > a big number k{e) £ N exists that 

He) 

,ii,i2,...,ik (-^J 

k=l ii,i2,---,ik 

It is clear that ^ oo as e — >■ +0. 

Let us take any e > and stop our process at the step k{e). I will represent all finite 
functions fk,h,i2,...,ik-ii'^') ^ k < k{e) as a difference of convex functions in the 
following way. 

Let us define the function [1] 

/fe,n,i2,...,ifc-i(^) = 1^8^ v€dfkiji2 ik i(^) ^•^*^'*i'*2.-,ifc-i(^) + (^'^ ~ ■^))' 

where z = {x,y) G ^^,z = {x,y) and dfk,i-i,i2,...,ik-i{^) subdifferential (see [1], [8] ) 

of the function fk,ii.i2 ifc-i(^) point z. The function fk,ii,i2,...,ik-ii'T') ^ convex 

function, because it is the maximum of the linear functions [1], [8]. 
It is clear that the difference 

/fc, 11,12, ~ /fc,ii,i2,...,ifc-i(^) ~ /fc,ii,i2,...,ifc-i(^) € -D 

is a convex function equal to zero on -Dj^ j2,...,ifc_i- Prom here I have 

fk,ii,i2,...,ik~i{^) ~ fk,ii,i2,...,ik-ii^) ~ fk,h,i2,---,ik-ii^) ^ ^^ 

The proof is completed. 

I will obtain the necessary and sufficient conditions when the sums 

kie) 

iz) (3) 

11,12, ■■■,ik-i 

and 

k(e) 

(z) (4) 

11,12, ■■■,ifc-i 

are bounded by a constant that does not depend on e,z = {x, y) G D. 
Let us take any r(-) G p{D). Let 

$(t) = /(r(t)) VtG[0,r,], 

where r(.) is written in the natural parametric representation. The parameter Tj. is the 
length of the loop r(.). 

I will prove that $(•) is Lipschitz with the constant L. In fact for any ti,t2 G [0,Tr] I 
have 

I Hh) - $(t2) 1=1 /(r(ti)) - f{r{t2)) \< L\\r{h) - r{t2)\\ <L\h-t2\. 
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Therefore [9] $(•) is almost everywhere differentiable in [0,Tj.]. I wiU define by Nr the set 
of the points where the function $(•) is differentiable in [OjTr]. 

I will prove that for any curve r(-) G p{D) there is a constant c{D) > such that 

V($';0,r,) <c(D). (5) 

Then the summations (3) and (4) are bounded on module by some constant that does 
not depend on e > 0. Since the function /(•) is continuous on D, it is sufficient to prove 
that one of the summations, (3) or (4), is bounded. 

The proof will be based on the following Lemma 2.1 and Lemma 2.2. 

Lemma 2.1 . The inequality 

yiQ';0,Tr) <ciiD,^) 

is true for any convex positively homogeneous functions with the degree 1 Q ^ ipiq) '■ 
— )■ 5ft and for any r(-) G p(-D) where Q(t) = ip(r{t)) for any t G [0,Tr],ci(D,xp) is a 
constant. 

Proof. Without loss of generality I will consider that 'ip{-) is a smooth function in 5ft^ {0} 
and {0} G -D. Let 

V'(r(t)) = max,e9^(o)(i;,r(t)) = {v{t),r{t)), v{t) G 9V^(0), 

where dtp(0) is the subdifferential of the function V'(-) at the initial point. As well I will 
consider that r(.) is a differentiable on t function. As soon as r(t) is the support vector of 
the set dil){G) at the point v{t), we have {v'{t),r{t)) = 0. 
I have 

I V'(r(ti))-V'(r-(t2)) 1=1 {v{tr)y{t^))-{v{t2)y{t2)) \ = \ {v{t,) - v{t2),r'{t-,))+ 

{v{t2),r'{t,)) - {v{t2),r'{t2)) \< \\v{t,)-v{t2)\\ \\r'{t,)\\ + \\r' (ti) - r' {t2)\\ \\v{t2\\ < 

\\v{ti) - V{t2)\\ + L{D) \ti-t2\. 

Here we took into consideration that ||r'(t)|| = 1. 
It follows from here that 

V(e';0,r^) < 2P{di;iO)) + L{D)Tr, 

where P(V'(0)) is the length of the curve bounding the convex set dtp{0) C and L(D) is 
the Lipschitz constant of the function tp{-)- Since all terms in the right side depend only 
on the set D and some constants, Lemma 1 is proved. □ 

Lemma 2.2 Let (x, y) — t- y) : 3f?^ — )■ 3^ 6e a continuous convex function and r(-) G 
p{D),t G [0, Tr]. Then there is a constant C2(D,/i) > such that 

V($i;0,r,)<C2(A/i) (6) 
where = /i(r(t)), t G [0, T^]. 



9 



Proof. Without loss of generality I will assume that is not negative on D, the 

initial point is the minimal point and also = (0, 0) belongs to the interior of the convex 
region in 3?^ with the boundary r(-). I will construct an homogeneous function -(/'(•) with 
the degree one, i.e. Tp{Xx,Xy) = X7p{x,y),X > 0, with values on the curve r(-). I 

will prove that ^{■) is the convex function. 
Let us consider the function 

fe{x,y) = fi{x,y) + e{\\ X ||2 + II y ||2),£ > 0. 

I will divide the segment [0, T^] into equal segments using division points {U}, i E 1 : J, 
and build the planes tTj in 5?^ containing the points (0,0,0), {r{ti), fg{r(ti))), 
(r(ii+i), /£(r(ti+i)), i G 1 : J. The parts of the planes vrj, i G 1 : J , are defined in their sec- 
tors (0, 0), (0, r{ti)), (0, r{ti+i) and built with the help of the vectors (0, 0, 0), {r{ti), fsiriU))), 
{r{ti-\-i), fe{r{ti+i))). They define a homogeneous many-sided function (■0£)j(r(-)) (refer 
[3], [7]). I will prove that all two-sided angles defined by the planes tTj, i € J, of the function 
are convex. 

Since any curve r(-) G p{D) can be approached by a smooth curve from p{D), without 
loss of generality we can assume that r(-) is a smooth curve with the derivative r'(-). 

Let us define the gradients of the planes tTj and tTj+i by VvTj and V-TTj+i respectively. 
I will use the theorem about midpoint according to which there is a point tm G 
such that 

dfe{r{tm))/dei = {Vnuei), 

where 

Cj = {r{ti+i) - r{ti))/ II r{ti+i) - r{ti) \\= r'{tm)- 
Similarly for the plane tTj+i and some point tc G [tj+i,fi+2] I have 

dfe{r{tc))/dei+i = (VTTj+i, Cj+i), 

where 

Ci+i = (r(tj+2) - r{ti+i))/ II r(tj+2) - r{ti+i) \\= r'{tc). 

The function f^{-) has a positive definite matrix of the second derivatives. Consequently, 
according to the well-known qualities of convex functions (see, [1] ) for sufficient big J and 
uniform subdivision of the curve r(-) by points ti we have 

dfe{r{tm))/dei < dfMtc))/dei+i, 

or 

(V7rj,ei) < (V7rj+i,ei+i). 

Let us note that the difference V-TTj+i — ViVi is perpendicular to the vector r(tj+i). From 
here and the inequality written above it follows that all two-sided angles tt, , vrj-i-i are 
convex. For J — t- oo we have the uniform convergence (^e)j(-) =^ {ip£){-) on D. 

As soon as the pointwise limit is equivalent to the uniform limit, then Vel') is the 
convex function. Also we have the uniform convergence ipei') =^ V'C') when e — > -|-0, i.e. 
^(•) is convex on D as well. 
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It is obvious that the gradients of the hnear functions with the graphs vrj , i G J , are 
bounded by a constant depending only on the region D and the function •) and 

V^(r(t)) = /i(r(t)) yte[0,Tr]. 

It is obvious that V'(") ") was built correctly with the help of the function •) and the 
chosen curve r(-). It follows from above that the function tp(-, •) is Lipschitz with a constant 
L(D, f). 
Let 

= Vie [o,r,]. 

Since 

v($;;0,r,) = v(^'i;0,r,), 

it follows from the Lemma 2.1 that 

V($i;0,r,)<C2(A/i). 

Hence the Lemma 2.2 is proved. □ 

It follows from the Lemma 2 that if /(•, •) is represented as a difference of two convex 
functions i.e. 

f{z) = Mz)-f2{z) \/zeD 

where •), i = 1, 2, are convex functions then the inequality (6) is true by necessity. In 
fact, I will introduce for any r(-) G piD) the following notations 

= fiir{t)),-^2{t) = f2{r{t)) Vi G [0,r,]. 

Since 

V($'; 0, Tr) < V($i; 0, T,) + y{^'^- 0, T,) 

the inequalities (5) and (6) are true by necessity [9]. 
I will prove, for instance, that the sum 

(t) 

11,12, 

is bounded by some constant on D independently from e where A;(e) — )• oo with e — >■ +0 
(in general case) when (5) is true for any curve r(-) G p(-D). 
I will introduce the functions for any r(-) G p(-D) 

4Mi,i2,...,i.-iW = /fc,n,i2,...,i,-i(Ki)) vtG[o,r,], vfc, 

I will prove the inequality 

kis) 

EvKn,i2,...,..-i;0,T.)<Ci (7) 
k=l 
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for any small e > O.Here Ci does not depend on e. 

In one-dimensional case the functions fk,ii,i2,...,ik-ii'^ ') ^'^^^ t)e linear in regions where 
they do not coincide with fk-i.ii,i2,...,ik-i ('■> ')■ Therefore the sum of variations of derivatives 
of the functions fk,ii,i2 ik-ii'^ ') ^° exceed the same variation of the function /(•)• 

I two-dimensional case the functions fk,ii,i2,...,ik-i{'->') be linear in regions 

where they do not coincide with fk-i,ii,i2,...,ik-i{': ")• I^or this reason the sum of variations 
of derivatives of the functions fk,ii,i2,---,ik-ii'i ') t)e greater than the same variation of 
the function /(•). I will prove that the total increase on k will be bounded above. 

The proof will be done from the contrary in the same way how boundcdncss of 
(3) was proved from (5). I will take segments in the regions -Dfc_i j2,...,ii_i) where 
/fc,ii,i2,...,H_i(-)-) do not coincide with fk-i,h,i2,...,ik-ii'^')^ slong which variations Vk of 
derivatives of fk,h,i2,...,ik-ii'i') maximal. But V]^ are not greater than variations of 
/(•) along some segments from -Dfc-i,ii,i2,...,ii-i- If the sum of Vk is not bounded above, 
then we can find out a curve or a sequence of curves from p{D) along which variations of 
derivatives of /(•) are not bounded above. It contradicts (5). Consequently, (7) is proved. 

I will suggest the contraposition. Namely, let the inequality (7) be true but the sum 
(3) is not bounded on D when e — > -|-0. 

I will choose for every function fk,h,i2,-,ik-ii' ^ ') ^ direction along that the variation 
of the derivative of the considered function fk,h,i2,...,ik-ii'j ') maximal. I will prove that 
the constant Ci not depending on e in (7) docs not exist for our suggestion. 

Let us denote by Ik,ii,i2,-,iki ^ segment of the region -DA;,ii,i2,...,ij_i along that the vari- 
ation of the derivatives of the function fk,h,i2,...,ik-ii'i is maximum among all segments 
of the region -Dfe,ii,i2,...,ii_i. 

Let us make the following procedure. Dividing the region D into convex compact 
subregions I will choose such a subregion where the function /(•,•) is not represented as 
a difference of two convex functions. Furthermore I will divide the chosen subregion into 
smaller subregions and chose one of these smaller subregions where the function /(•,•) is 
not RDC function. I continue this process until I obtain a point M such that the function 
/(•, •) is not RDC in a neighborhood of the point M. Further I will divide the set D into 
sectors that have a common point M, and choose a sector, where the function /(•, •) is not 
RDC function. The set of the chosen sectors shrinks to a ray I with the summit at the 
point M. It is obvious that the function /(•, •) is not RDC in any cone K with the summit 
at the point M if I G intK. 

From the procedure described above I will produce one of the following two situations: 

a) the sum (3) has the infinite variation along the direction I when e — >■ -|-0; 

b) the sum (3) has the infinite variation when s — >■ -|-0 along the direction rj that is 
perpendicular to the direction I. 

The ) and b) mean that the infinite sum (3) of the convex functions is an 

unbounded convex function in the directions Z or r/ correspondently. 

Consider the case a). Take any cone K with non-empty interior and I G intK. I will 
consider the functions fk,ii.i2,...,ik-i('^ ') with the supports D^ i^ j,, from the cone K. 

Let the slope angles between the segments Ik,h,i2,...,ik-i ^^'^ the direction / for large k 
do not exceed 7r/2 — (5 for some (5 > 0. Then I will choose a subgroup of the segments 
Ik,ii,i2,...,ik-i ttiEi't can be intersected by a curve r(-) G p{D) that has the slopes at the 
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points of intersections with the segments Ik,ii,i2,---,ik-i exceeding 7r/2 — 6,6 > 0. Since 
the cone K is arbitrary and contains the ray I I can assume that r'{t) — > I when t ^ Tr. 
As soon as all functions fk,ii,i2,...,ik-i{'^') Lipschitz uniformly, their variations of 

the derivatives along a segments that is close to the segment Ik,ii,i2,...,ik-i ^^^^ 
close to each other as well. The curve r(.) will consist of the segments Sk in the regions 

Let a curve r(-) G p{D) exist along that the variation of the derivative of the sum (3) 
increases to infinity when e — t- +0 (refer to the fig. 3). (It is possible the case when the 
regions Dk,i^,i2,...,ik-i include each other but it does not change the discussion). I obtain 
the contradiction with the inequality (5). 

The case is possible when I have several groups of the segments {Ik,i\,i2,...,ik-i\i 
for each of these groups a curve ri(-) G p{D) exists that 

V«^;0,T,,) = Ci, r[{t)^l 

when i — )• T^., i is the natural parameter, and 

i 

where ^ni') is the value of the infinite sum (3) for e — t- +0 on the curve ri(-). 

I will build a curve r(-) G p(-D) for this case in the following way. 

The curve r(-) has to consist of enough number of segments from each group (or close 
to them) so that 

^{KvUi^Ui+i) = Ci-£i 

where is the value of the parameter t for the i-th group of the segments for the natural 
parametric representation of the curve ri{-) , Sj is a positive number so that 

^ £i = £■ 

i 

It is obvious that such a curve r(-) can be always built. It will consist from the curves 
ri(-). For that it is necessary to do smooth passage from one curve ri(-) to another r(-) 
without going out of the set p{D). Since r^(t) — >■ Z for t — >■ T,. and all i, this procedure can 
be always done. But then 

V«;0,r,) > ^ V($;^;0,T,J = ^(q - e,) = ^ - ^ = go. 

i i i i 

where 4>r(0 is a value of the infinite sum (3) for e — >■ +0 on the curve r(-). 

As soon as, according to the algorithm, the sum of the variations (7) along any segment 
of the curve r(-) G p{D) is equal to the variation of the initial function along the same 
segment, then we have 

V($';0,T^) = oo 

where $(•) is calculated at the constructed curve r(-) that contradicts (5). 
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If the angles of the slopes of the segments Ik,ii,i2,-,ik-i with the ray / go to 7r/2 when 
A; — >■ oo , then the variation of the sum (3) is infinite in the direction rj when e — >■ +0 (the 
case b). 

Let us consider the case b). 

We will construct the curve r^(.) from the segments Sk in the regions I?fc,ji,i2,...,ife_i- 
All the segments Ik,h,i2,-,ik-i divided into such groups {m} of segments that 

can be intersected by a curve r^(.) G p{D) for that 

where is a parameter of the curve rm{-) for the natural parametrization at the point 
M. The curvatures of the curves rm{-) tend to infinity when r — >■ Tr^. Obviously, it can be 
always done by dividing the set of the segments Ik,h,i2,...,ik-i subsets with required 

qualities (refer to the fig. 4). 

Moreover, the angles a^^, the curves rm{-) and the groups of the segments {Ik,ii,i2,-,iki}m 
can be chosen in such a way that 

^TooV(^^">5 0,r,^) = oo. (8) 

The construction of the curves rm(-) G p{D), for that (8) is true, is fulfilled in the same 
way as in the case a). For that it is necessary to intersect by the curve r^(-) G p(-D) 
enough number of segments from each group that 

and 

lim Cm = oo, 

m—^oo 

where [trm i '^vm+i] the segment of values of t of the natural parametric representation for 
the m -th group of the segments. We can always construct such a curve rm{-) G p{D). 
The inequality written above will be satisfied as soon as 

1. Increasing the curvatures of the curves at the point M the curves rm{-) cross an 
increasing number of segments Ik,n,i2,...,iu^ for m ^ oo with acute angles am > 0. 

2. A number of segments is infinite in any cone K with the vertex M and the inner 
vector 1. 



3. The sum of projections of their lengths in the direction -q is infinite. 

The described procedure gives the finite functions /fc,n,j2,.--,«fc-i("' ")• The sum of the vari- 
ations of the derivatives of these functions along the segments Sk of the curve r(.) is equal 
to the variation of derivatives of the function $(t) = f{r{t)) along the same segments. As 
a result we have 

lim V($';0,T,^) = oo, 
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where $(.) is calculated along the constructed curves rjyi(.). We come again to the contra- 
diction with (5). 

Two operations described above in this paper allow us to come to finite or countable 
number of the finite functions from 9 as long as the maximums of the modules of the 
functions fk,ii,i2,...,ik-i{'i ') 8° to zero in D when k — > oo. Since the sum of the variations 
of the derivatives of the functions fk,ii,i2,...,ik-ii': ') along any segment from D is bounded 
above by the constant c{D, /) , Theorem 1 is proved. □ 

Remark 2.1 The arguments about the chosen angles and the curves rm{-) are anal- 
ogous to the following. 

Let us have the divergent series 



Here ai is the analogue of the variation of derivatives along a segment I k.i^. 12, ■■■,ik-i 

is the analogue of the cosine of an angle between a curve ri and a segment Ik,h,i2,-,ik-i 

the point of intersections. 

3 Geometric interpretation of the theorem 1 

I will give another geometrical interpretation of the Theorem 1. Let me introduce a turn 
of the curve r(-) on the graph T j = {{x,y,z) £ \ z = f{x,y)}. 

Consider on Tf the curve R{t) = {r{t), f{r{t))) where r(-) G p{D). As the function 
/(■,■) is Lipschitz, the derivative R'{-) exists almost everywhere on [0, T^], which I will 
denote by t(-) = R'{-). 

Definition 3.1 The turn of the curve R{-) on the manifold Tf is called the value 



Thus the turn Or of the curve R{-) is the supremum of the sum of all angles between the 
tangents r(t) for t G [OjT^]. It is obvious that the value Or is equal to the integral 




We can always choose a decreasing sequence {A}, A — >-i-).oo 0, that 



m 




1=1 





for any fiat curve r(-) represented in the natural parametric way where k{s) is the curvature 
of the curve r(-) at the point s G [0, T^] i.e. in this case this definition coincides with the 
usual definition of the turn of the curve [10] . 
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Theorem 3.1 . For any Lipschitz function z f{z) : D ^ R to be RDC function on 
the compact set D E it is necessary and sufficient that a constant 02(0, f) > exists 
for all r(-) G p(-D) such that the turn of the curve R{-) on the Tf was bounded by the 
constant C2{D,f) > depending only on D (and the function f) i.e. 

Or<C2{D,f) WepiD). (9) 

Proof.Necessity. Let /(•,•) be a RDC function. I will prove that the inequality (9) is 
true in this case. Let us use the inequality following from the triangle inequality 

\\T{ti)/\\T{ti)\\-Tit,.,)/\\Tit,.4\\<\^^^ 

I f'MU))/^! + f?{rm - //(r(t,-i))/Vl + /f I - 

Since 1 < y^l + fp{r{ti)) < vTTX^ for all ti G [0, T^] , it is obvious that the number 
C3 > 1 exists, for that the inequality 

\\r'{t.)/^l + f[Hr{t^) - r'(t,_i)/^l + //2(r(t,_i))|| < c^\\r'it,) - r'{t,.,)\\ (10) 
is true. The inequality 

I fMU))/yJl + /f (r(t,))-//(r(t,_i))/^l + /f (r(t,_i)) |<| //(r(t,))-//(r(t,_i)) | (11) 
follows from properties of the function 6{x) = x/VY+lc^. We have from (10) and (11) 
sup Yl l|T(ii)/||r(ii)||-T(ii_i)/||r(t,_i)|||| <C3(V(||r'||;0,r,) + V($';0,T,)). (12) 

{ti}eNr i 

Since the function /(•,•) is RDC, the inequality 

V($';0,T,) < ciiD) 

is true according to Theorem 1. The inequality (9) follows from the inequality shown 
above and from (12). The necessity is proved. 

Sufficiency. Let the inequality (9) be true. I will prove that /(•,•) is the RDC 
function. I will use the inequality 

||r(t,)/||r(tOI|-r(t^-i)/||r(t,_i)|||| >| fl{r{ti))/yjl + fl{r{ti))-fl{r{ti.,))/^l + flHr{U-i)). 

(13) 

From the qualities of the function 9{x) = x/\/l + and from the inequality ||/'(-z)|| < L 
for z G D it follows the existence of a constant C4(L) > 0, for that 

I f'Mti))/^l + fp{r{ti)) - //(r(t,_i))/^l + /f (r(t,_i)) > c, \ f'Mt,)) - fi{r{ti^i)) \ . 
Prom here and (13) we have 

C2> sup ^||r(ti)/||T(t»)||-r(t,_i)/||r(t,_i)|||| >C4V($';0,r,). 

{ti}(ZNr i 

It follows from Theorem 1.1 that /(•) is a RDC function. Sufficiency is proved. □ 
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